Periodically rippled graphene: growth and spatially resolved electronic structure 
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We studied the growth of an epitaxial graphene monolayer on Ru(0001). The graphene monolayer 
covers uniformly the Ru substrate over lateral distances larger than several microns reproducing the 
structural defects of the Ru substrate. The graphene is rippled with a periodicity dictated by the 
difference in lattice parameter between C and Ru. The theoretical model predict inhomogeneities 
in the electronic structure. This is confirmed by measurements in real space by means of scanning 
tunnelling spectroscopy. We observe electron pockets at the higher parts of the ripples. 

PACS numbers: 73.20.-r, 68.37.Ef, 68.55.-a, 81.05.Uw 



The possibility to produce single layers of graphene 
0, [2] has opened a fascinating new world of physical 
phenomena in two dimensions. Graphene has already 
shown that its charge carriers are massless Dirac fermions 
[!, EL HI El and that it displays an anomalous integer 
Quantum Hall Effect 0, 3, Q even at room temperature 
[8[ . Systems made up of a few graphene layers can also 
be grown on a SiC substrate [9|, ll0[. Recently it has 
been shown that free standing isolated graphene layer is 
intrinsically corrugated [ll|. These modulations in the 
heights of the graphene layers may be related to the 
charge inhmogeneities observed in nominally undoped 
samples [ 12] (see[13j). 

Ultra-thin epitaxial films of graphite have been grown 
on solid surfaces for quite some time [14J. Even the 
growth of " monolayer-graphite" films onto several sub- 
strates by Chemical Vapor Deposition have been reported 
some time ago, but the degree of characterization of the 
films was hampered by the existing experimental limita- 
tions fist - 

In this letter we report on a method to fabricate highly 
perfect, on the scale of microns, periodically rippled 
graphene monolayer and islands on Ru(0001) under Ultra 
High Vacuum (UHV) conditions. We characterized by 
means of scanning tunnelling microscopy /spectroscopy 
(STM/STS) the perfection at atomic scale and the local 
electronic structure of the periodically rippled graphene 
monolayer. The periodicity of the ripples is dictated by 
the difference in lattice parameters of graphene and sub- 
strate, and, thus, it is adjustable. We observe inhomo- 
geneities in the charge distribution, i.e electron pockets 
at the higher parts of the ripples. This inhomogeneity in 
the electronic structure can be understood with the help 
of a tight-binding model. The potential associated with 
the rippled structure induce a charge transfer from con- 
duction to valence bands for some atoms and the opposite 
in the others. Finally we studied graphene nanoislands 
that display hexagonal shape with atomically resolved 
zigzag edges, whose characteristic edge states allows us 



to determine the doping of the graphene layer. 

The experiments have been carried out in a UHV 
chamber with base pressure of 4xl0 -11 Torr that con- 
tains a variable temperature STM and a rear view 
Low Energy Electron Diffraction (LEED) optics. The 
Ru(0001) crystal was cleaned by cycles of Ar + sputtering 
and annealing followed to oxygen exposure and heating 
to high temperature [16[ . The graphene layers were pro- 
duced either by controlled segregation of C from the bulk 
of the substrate or by thermal decomposition at 1000 K 
of ethylene molecules pre-adsorbed at 300 K on the sam- 
ple surface. The W tips were routinely cleaned by ion 
bombardment and annealing. The dl/dV curves were 
obtained by numerical differentiation of the I(V) curves. 
We model the graphene electronic bands by a nearest 
neighbor tight-binding model with parameter t=3eV, a 
periodic shift of the carbon p z levels, and a finite broad- 
ening due to the hybridization with the metallic sub- 
strate. In order to calculate the density of states, we 
use a (30x30) unit cell, and a sum over 6 special points 
in the irreducible sector of the Brillouin zone. Hence, 
the total number of states included in the calculation is 
10800. This is sufficient to resolve local changes in the 
electron density in the order of 10 -2 -10 -3 carriers per 
carbon atom. 

The epitaxial layer of graphene covers completely the 
surface of the single-crystal Ru substrate over distances 
larger than a micron (see [17]) and presents a triangular 
periodicity of 2.4 nm (Fig. 1(a)) that is due to the co- 
incidence lattice of graphene and Ru, i.e. the lattice of 
graphene is incommensurate with the underlying lattice 
of Ru, with a size relation that implies that 10 carbon 
honeycombs (0.246 nm) will adjust almost exactly with 
9 Ru-Ru interatomic distances (0.27 nm). The weakly 
interacting (see below), laterally undistorted graphene 
structure rides on top of the lattice of the substrate, re- 
sulting in some C atoms being somewhat higher than 
others. The graphene layer, thus, is rippled in a peri- 
odic fashion, with a periodicity dictated by the difference 
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in lattice parameters. Variations of this approach using 
other experimental conditions can be used to produce 
graphene monolayers on other single-crystal substrates 

With the STM images we characterize the type of de- 
fects that appear in the epitaxial graphene layer. The 
monoatomic steps and dislocations of the substrate are 
reproduced, but there are also upper part of the ripples 
that are weaker or even missing completely in the STM 
images, as can be seen in Fig. 1 (a) and in ref. [25]. Fig. 
1 (b) shows a high resolution topographic STM image of 
the epitaxial graphene layer, revealing its atomic struc- 
ture. Unlike STM images of graphite [19], which show 
normally only one of the two C atoms in the surface unit 
cell, the honeycomb structure of graphene is clearly re- 
solved in the upper part of the ripples with its 0.14 nm 
C-C distance. The Fourier transform of larger atomi- 
cally resolved STM images shown as inset in Fig. 1(b) 
indicates that the ripples constitute a (10x10) superlat- 
tice with respect to the C lattice. The (9x9) coincidence 
lattice with the Ru(0001) substrate is revealed in the cor- 
responding LEED pattern (not shown). 

The apparent vertical corrugation of the rippled 
graphene monolayer, as seen with STM, changes with 
the tunnelling voltage, ranging from 0.1 nm (at a sample 
voltage of -1 V) to 0.02 nm (at +1 V). This indicates 
the importance of electronic effects. The C-C apparent 
atomic corrugation inside the honeycomb unit cell is only 
of the order of 0.005 nm (Fig 1(d)). The actual geomet- 
ric corrugation of the rippled layer can be estimated to 
be below 0.02 nm from the observed lack of diffracted 
peaks for a beam of thermal He atoms scattering off the 
graphene surface [2o| . 

The bonding with the substrate occurs through the 
hybridization of the C 7r-states with the Ru d states. 
Photoelectron spectroscopy shows clearly that the layer 
bonding is not carbidic. Angular Resolved Photo Elec- 
tron Spectroscopy data show that the graphene bands in 
graphene/Ru(0001) are similar to the ones of graphite, 
but rigidly shifted down in energy [2l|. The bottom of 
the 7r band at the center of the Brillouin Zone is shifted 



down by 1.8 eV with respect to graphite [2l|. The small 
energy shift of the Cls core level with respect to graphite 
indicates that the charge transfer from the substrate is 
small, but not negligible, i.e. the graphene layer is doped 
with electrons by the substrate. 

The graphene-Ru interface is atomically well defined. 
The average interlayer distance for the closely related 
incommensurate graphene/Pt(lll) system has been de- 
termined to be rather large (0.37 nm) by means of ten- 
sor LEED analysis [l5j]. First principles calculations in 
the similar graphene/Ir(lll) system have indicated that 
the average binding (0.2 eV/C atom) is barely strong 
enough to correspond to chemical bond formation [22| . 
The incommensurate nature of the graphene monolayer is 
an additional indication that the interaction of graphene 
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FIG. 1: (Color on line) : (a) 76 nm x76 nm STM image of 
graphene/Ru(0001) showing the decoration of a screw dislo- 
cation and a monoatomic step from the substrate. There are 
also some defects on the rippled structure, (b) 6.5 nmx6.5 
nm atomically resolved image of graphene/Ru(0001). The 
image was taken with a sample bias voltage of V s =-4.5 mV 
and a tunnel current of It =3 nA. The image is differentiated 
along the X direction in order to see the weak atomic cor- 
rugation superimposed to the ripples. The inset reproduces 
the Fourier transform of the image showing the (10x10) pe- 
riodicity of the rippled graphene layer. The larger hexago- 
nal pattern corresponds to the C-C distances and the smaller 
spots to the periodic ripples, (c) Corresponding structural 
model, (d) Line profile marked with an arrow in panel (b). 
The atomic corrugation is around 5 pm. 



with the substrate is rather weak. 

Fig. 2(a) shows spatially-resolved dl/dV tunnelling 
spectra, which are roughly proportional to the Local Den- 
sity of States (LDOS), recorded on top of the "high" and 
"low" regions of the corrugated graphene layer. Panel (b) 
shows the corresponding calculation. The experimental 
tunnelling spectra recorded at different spatial positions 
are obviously different: the occupied LDOS is system- 
atically larger in the "high" areas of the rippled layer, 
while the empty LDOS is larger in the "low" parts. The 
differences are robust enough to survive at 300 K. There 
are also weaker features at both sides of the Fermi level, 
separated from each other by 0.3 eV and located at about 
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FIG. 2: (Color on line) (a) dl/dV curves recorded at 300 K 
on the higher (black curve) and lower (red curve) areas of 
the rippled graphene monolayer. The inset shows the topo- 
graphic image measured simultaneously, (b) Calculations of 
the spatially resolved LDOS for an isolated, (10x10) peri- 
odically corrugated graphene layer. The spectra correspond 
to the high (black curve) and low (red curve) regions of the 
rippled graphene layer. 



the same energies in both regions of the ripples. 

The electronic structure has been simulated by taking 
into account that the rippled graphene layer is not too 
distorted. The model calculations have been performed 
for an isolated graphene layer in which the effect of the 
substrate has been considered to result in: i) a shift of the 
Dirac point by -0.3 eV due to doping; ii) the introduction 
of a finite lifetime caused by hybrization of the ir orbitals 
with the band of the substrate; and iii) a (10x10) peri- 
odic potential that changes between -3V and 3V/2,where 
V=-0.3eV to account for the periodic structural ripples. 
This value for the potential gives the better agreement 
with the experiments. For higher values of the potential 
the asymmetry in the calculated DOS is bigger than the 
one found in the experiments. The ir band in graphene 
has a total width of W~6t, where t~3 eV is the hopping 
between tt orbitals at nearest neighbor C atoms [23[. 

In agreement with the experiments, the calculations 
show that the occupied LDOS is larger on the "high" 



regions of the superlattice, where the potential is at a 
minimum, while the empty LDOS is larger at the "low" 
regions of the graphene layer (see lower panel in Fig. 3). 
This effect is extremely robust and indicates that the va- 
lence band is depleted in the low portions of the ripples, 
while the conduction band is depleted at the high parts 
of the ondulated graphene layer. The calculations also 
reproduce the weaker features at both sides of the Fermi 
energy that reflect the existence of the periodic super- 
lattice, which induces a folding of the graphene bands 
in the new (10x10) Brillouin Zone. The reconstruction 
consider here does not superimpose the two inequivalent 
corners of the Brillouin zone and it does not open a gap at 
the Dirac energy [24| . The spectrum is split on subbands 
separated by gaps, away from the K and K' points. 

The periodic charge inhomogeneities in the graphene 
layer can be visualized directly in the real space by imag- 
ing the spatial distribution of dl/dV close to the Fermi 
energy. Fig. 3 shows the spatial distribution of the 
LDOS below and above the common Fermi level (which 
corresponds to the shifted neutrality level of the doped 
graphene layer). The experimental images are in the up- 
per row at the left and the right of the corresponding to- 
pographic image. The bright regions correspond to larger 
LDOS. For a complete set of spatially resolved dl/dV 
maps versus sample bias voltage see ref. [251 ]. 

These inhomogeneities in the charge distribution pro- 
duced by the periodic ripples are probably also present 
in the structural corrugations inherent to free standing, 
isolated graphene monolayers [Tlj . 

The growth method described here can also be used 
to produce different graphene nanostructures, such as 
nanowires on vicinal surfaces or nanometer wide islands 
as illustrated in Fig. 4(a). The apparent step height of 
the island is 0.15 nm (from the Ru surface to the lower 
part of the ripples) showing that the graphene layer is 
indeed only one monolayer high. The inset shows that 
the islands are truncated hexagons with straight edges 
of a single structural type. It is remarkable that the 
periodically rippled coincidence lattice goes right to the 
steps of the islands. In fact in some cases the island step 
cuts the (10x10) arrangement. This confirms that the 
overlayer of graphene is indeed weakly coupled to the Ru 
substrate. Further zooming into the island reveals that 
the geometry of the steps is of the zigzag type. 

The relative orientation of the graphene layer with re- 
spect to the underlying Ru lattice can be determined by 
resolving simultaneously both atomic lattices. In order 
to facilitate the resolution of the weakly corrugated Ru 
lattice, CO has been adsorbed at saturation in the parts 
of the Ru surface not covered by graphene. The geome- 
try of the resulting (\/3 x y/3 R30°) structure of CO on 
Ru(0001) is well known [26[, and, thus, the relative ori- 
entation of graphene and Ru lattices represented in Fig. 
1 has been obtained. 

Fig. 4(b) shows local tunnelling spectra recorded when 




FIG. 3: (Color on line) The left and right images in the up- 
per panel are maps of dl/dV at -100 meV and + 200 meV 
and reflects the spatial distribution of the LDOS below and 
above the Fermi level, respectively, for an extended graphene 
layer on Ru(0001). The central image shows the topographic 
image recorded simultaneously. The lower panel shows the 
corresponding calculations for the spatially resolved LDOS 
for a (10x10) periodically corrugated graphene layer. The 
arrows indicate the (10x10) unit cell. 



FIG. 4: (Color on line) (a) 10 nmxlO nm atomically resolved 
STM image of a graphene island on Ru(0001). The image was 
recorded at V s =-4.5 meV and It =3 nA. The inset shows an 
image of the whole island with a lateral size of 47 nm. The 
images are differentiated in X direction, (b) Spatially resolved 
tunnelling spectra measured on the high (black curve) and 
low (red curve) areas of the ripples close to the edge of the 
island and on clean Ruthenium (blue curve) . The Ruthenium 
spectra was offset for clarity 



the surface is only partially covered with graphene. The 
advantage of this situation is that one can take spectra in 
the clean Ru patches of the surface, where only the sur- 
face state of Ru should be detected, and, thus, spurious 
features from the tip DOS can easily be detected. The 
spectra recorded at the center of the islands are undis- 
tinguishable from the one obtained on extended graphene 
layers (see Fig. 2). The spectra taken at the edges, how- 
ever, show a new peak at -0.4 eV, particularly on the 
"high" areas of the ripples, which might be related to the 
theoretically predicted edge states [27[. The spectrum 
recorded on the clean Ru patches show only the surface 
state around the Fermi energy observed on Ru(0001) pH ]. 

In summary we have grown graphene monolayers and 
islands on Ru(0001). The graphene monolayer presents 
a periodically rippled surface due to the difference in lat- 
tice parameter between graphene and the metallic sub- 
strate. The periodic ripples produce a spatial charge re- 



distribution in the graphene. This have been measured 
with spatially resolved dl/dV maps and confirmed with 
a theoretical model. The new periodicity also induce the 
opening of a series of minigaps by the additional periodic 
potential. The presence of these minigaps is expected to 
give rise to new phenomena at low temperatures in the 
presence of high magnetic fields. We atomically resolve 
the zig-zag edges of graphene islands and we use the po- 
sition in energy of the edge state to determine the doping 
of the graphene layer. The methods described here can 
be implemented on many other single crystal substrates, 
giving rise to a series of graphene monolayers with dif- 
ferent, substrate-dependent, periodic corrugations and, 
thus, opening the possibility to systematically test the 
electronic properties of controlled, charge inhomogeneous 
graphene layers. Methods to transfer the present finding 
to experimental systems adequate for lateral transport 
measurements are in progress. 
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